Extra- galactic magnetic fields and the second knee in the cosmic-ray spectrum 
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Recent work suggests that the cosmic ray spectrum may be dominated by Galactic sources up 
to ~ I0 17 ' 5 eV, and by an extra-Galactic component beyond, provided this latter cuts off below the 
transition energy. Here it is shown that this cut-off could be interpreted in this framework as a 
signature of extra-galactic magnetic fields with equivalent average strength B and coherence length 
l c such that By/ki ~ 2 — 3 ■ 10~ 10 G-Mpc 1,/2 , assuming l c < r% (Larmor radius at < 10 17 eV) and 
continuously emitting sources with density 10 -5 Mpc -3 . The extra-Galactic flux is suppressed below 
~ f 17 eV as the diffusive propagation time from the source to the detector becomes larger than the 
age of the Universe. 

PACS numbers: 98.70.Sa, 98.65.Dx 



I. INTRODUCTION 

Recent developments, both experimental and theoret- 
ical, have significantly broadened the landscape of ultra- 
high energy cosmic ray phenomenology. The High Res- 
olution Fly's Eye experiment has reported the detection 
of a high energy cut-off ~ 10 20 eV 0, as would be ex- 
pected from a cosmological population of sources. This 
experiment has also observed that the chemical compo- 
sition is dominated by protons down to ~ 10 18 eV, and 
by heavy nuclei further below, in agreement with prelimi- 
nary KASCADE data Q • This and the steepening of the 
cosmic-ray spectrum at the "second knee" ~ 10 17 5 eV 
suggest the disappearance of the low-energy (heavy nu- 
clei) component and the nearly simultaneous emergence 
of a high-energy (proton) component. On theoretical 
grounds, it has been observed by Berezinsky et al. Q that 
a cosmological distribution of sources producing a single 
powerlaw could fit the high energy part of the cosmic- 
ray spectrum from the second knee up to the cut-off at 
10 20 eV, including the dip of the "ankle" - 10 18 5 eV. In 
light of these results, it is thus tempting to think that 
the cosmic-ray spectrum consists of only two main com- 
ponents: one Galactic, another extra-galactic, with the 
transition around the second knee. 

There are alternative views, admittedly, in which the 
Galactic component dominates the all-particle spectrum 
up the ankle [4j ; this latter feature would then mark the 
emergence of the extra-galactic component rather than 
the signature of pair production as in Ref. This is- 
sue will be hopefully settled by ongoing and future cosmic 
ray experiments, through more accurate composition and 
anisotropy measurements. The discussion that follows as- 
sumes that the interpretation of Berezinsky and cowork- 
ers Q is correct, namely, that the transition between the 
Galactic and extra-galactic cosmic ray component arises 
at the second knee. 

This model then requires to impose a low-energy cut- 
off on the extra-galactic spectrum around 10 18 eV in or- 
der to not overproduce the flux close to the second knee. 



The exact position of this cut-off as well as the spectral 
slope below it must be tuned to how the Galactic com- 
ponent extends above the knee p|. 

The objective of the present work is to demonstrate 
that this cut-off could be interpreted as a signature 
of extra-galactic magnetic fields with average strength 
B and coherence length l c such that B\fQ ~ 2 — 3 • 
10 -10 G-Mpc 1 / 2 , assuming l c < Larmor radius at 
< 10 17 eV and continuously emitting sources with den- 
sity 10 -5 Mpc -3 . In this picture, the extra-Galactic spec- 
trum shuts off below 10 17 eV as the diffusion time from 
the closest sources becomes larger than the age of the 
Universe. The first knee is viewed here as the maximal 
injection energy for protons at the (Galactic) source. 

The existence of extra-galactic magnetic fields is of im- 
portance to various fields of astrophysics, including ultra- 
high energy cosmic ray phenomenology, but very little 
is known on their origin, on their spatial configuration 
and on their amplitude The upper limits on B\/T c 
from Faraday observations lie some two orders of mag- 
nitude above the value suggested here. In the present 
framework, experiments such as KASCADE-Grande @ 
could probe these magnetic fields thanks to accurate mea- 
surements of the spectrum and composition in the range 
10 16 -► 10 18 eV. 



II. PARTICLE PROPAGATION 

The main effect of extra-galactic magnetic fields on 
~ 10 17 eV particles is as follows. In a Hubble time, 
these particles travel by diffusing on magnetic inho- 



mogeneities a linear distance d 
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65 Mpc (Zscatt/1 Mpc) 1 / 2 , with i SC att the scattering length 
of the particle. If d is much smaller than the typi- 
cal source distance, the particle cannot reach the de- 
tector in a Hubble time; since i sca tt> hence d, increases 
with increasing energy, this produces a low-energy cut- 
off in the propagated spectrum. Current data at the 
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highest energies, notably the clustering seen by vari- 
ous experiments, suggests a cosmic-ray source density 
n ~ 10~ 6 — Mpc -3 0, which corresponds to a source 
distance scale ~ 50 — 100 Mpc. Hence Z sca tt 0.3 — 1 Mpc 
at 10 17 eV would shut off the spectrum below this en- 
ergy @. 

To be more quantitative one has to calculate the propa- 
gated spectrum and compare it to the observed data. The 
low energy part (< 10 18 eV in what follows) of the extra- 
galactic proton spectrum diffuses in the extra-galactic 
magnetic field since the scattering length Z S catt ^ d (d 
typical source distance). In contrast, particles of higher 
energies (> 10 19 eV in what follows) travel in a quasi- 
rectilinear fashion, meaning that the total deflection an- 
gle 9 lms < 1, since Z Bcatt > l loss , where Zi oss < 1 Gpc at 
E > 10 19 eV is the energy loss length (which gives an up- 
per bound to the linear distance across which particles 
can travel). 

In the diffusion approximation, the propagated differ- 
ential spectrum reads (see the Appendix): 
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The sum carries over the discrete source distribution; ri 
is the comoving distance to source i. Note that a factor 
(oo/o e ) -3 in Eq. (|A.2(I . with ao and a D the scale fac- 
tor at observation and emission respectively, has been 
absorbed in defining a comoving source density; the re- 
maining factor (ao/a e ) « dE g /dE, see below. The func- 
tion E g (t, E) defines the energy of the particle at time t, 
assuming it has energy E at time to. This function and 
its derivative dE g /dE can be reconstructed by integrat- 
ing the energy losses 0- In the Appendix, it is shown 
that Eq. (JIJ provides a solution to the diffusion equation 
in the expanding space-time under the assumption that 
the energy loss of the particle is dominated by expan- 
sion losses, which is found to be a good approximation 
for particles with observed energies E < 10 18 eV. In this 
case, dE g /dE « E g /E w ao/a c . Although photopair 
and photopion production losses on diffuse backgrounds 
are negligible with respect to expansion losses for most of 
the particle history, they set the maximum linear distance 
(hence the maximum time) across which a particle can 
travel. The time integral in Eq.JTJ is indeed bounded by 
the maximal lookback time t at which E g (i, E) — £" max 
and by the minimal lookback time to — t~ Z sca tt/c neces- 
sary to enter the diffusing regime, taken to be the solution 
of r(t) = X(t, E), where r(t) = f*° dt'/a(t') is the comov- 
ing light cone distance. The (comoving) path length A is 
defined in Eq. (fO)l as A 2 = f t ° dto -1 (t)D [a e E e /a{rj)], 
with E e = E g (t e ,E) the energy at injection. The physi- 
cal meaning of A is that of a typical distance traveled by 
diffusion, accounting for energy losses. 

The injection spectrum extends from some minimum 
energy (< 10 16 eV in the present model) up to i? max = 
10 22 eV (the exact value is of little importance here). 
The function Q(E g ) = K(E g / -E max )~ 7 gives the emission 



rate per source at energy E g , K a normalization factor 
such that J dE EQ(E) = L, with L the total luminosity, 
which is assumed to scale as the cosmic star formation 
rate from This theoretical star formation rate his- 
tory agrees with existing data at moderate redshifts and 
provides an argued prediction for higher redshifts. It also 
fits in nicely the constraints of the diffuse supernova neu- 
trino background |10| which would be violated by more 
steeply evolving star formation rates. The choice of the 
star formation rate is not crucial to the present analysis 
since the exponential cut-off due to the magnetic horizon 
dominates the effect of the star formation history on the 
low energy part of the spectrum. 

Only continuously emitting sources are considered 
here, although the effect of a finite activity timcscale for 
the source is discussed further below. The cosmologi- 
cal evolution of the magnetic field has been neglected for 
simplicity; if the diffusion coefficient depends explicitly 
on time t the solution Eq. remains exact. 

At higher energies, the propagated spectrum is given 
by: 
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and tj in Eq. is related to r.; by r, = J t " dt'/a(t'); r.; 
should not exceed X(ti,E), beyond which point motion 
must have become diffusive. 

The Galactic cosmic ray component is modeled as fol- 
lows. Supernovae are accepted as standard acceleration 
sites, yet it is notoriously difficult to explain accelera- 
tion up to a maximal energy ~ 10 18 eV. Thus it is as- 
sumed that the knee sets the maximal acceleration en- 
ergy for Galactic cosmic rays: in this conservative model, 
the spectrum of species i with charge Z takes the form 
j z (E) « (E/E z )-^ exp(-E/E z ), with 7l - 2.4 - 2.7 a 
species dependent spectral index, Ez = Z x E p the loca- 
tion of the knee, with E p m 2 ■ 10 15 eV 2j. This scenario 
is consistent with preliminary KASCADE data. 

Datasets from the most recent experiments are consid- 
ered here: KASCADE 10 15 -> 10 17 eV @, Akeno 10 15 -> 
10 186 eV [ll|, AGASA 10 18 5 -> 10 20 5 eV [12 , HiRes 
1017.3 _^ KPeV Q and Fly's Eye 10 17 3 -> 10 20 eV [l]j|. 
Akeno is not recent but it is the only experiment whose 
data bridge the gap between the knee and the ankle. 
These experiments use different techniques and their re- 
sults span different energy ranges, hence the data do not 
always match. A clear example is the discrepancy be- 
tween HiRes and AGASA at the highest energies. In 
the following, the high energy datasets have been split 
in two groups: one with HiRes and Fly's Eye, the other 
with Akeno and AGASA. The flux of the extra-galactic 
component was given two possible normalization values 
in order to accomodate either of these datasets, while the 
flux of the low energy part is scaled to the recent KAS- 
CADE data. A more robust comparison of this model 
with the data will be possible with the upcoming results 
of the KASCADE-Grande experiment covering the range 
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Energy [eV] 

FIG. 1: Open squares: Akeno; filled circles: KASCADE; 
filled diamonds: AGASA; filled squares: HiRes-2; open tri- 
angles: Hires-1; open diamonds: Fly's Eye. Fit of the Galac- 
tic (low-energy dot-dash line) and extra-galactic (high-energy 
dot-dash) to cosmic ray data. Total flux: solid line; dotted 
lines: upper 75 th and lower 25 th percentiles for the prediction 
of the extra-galactic flux. 

~ 10 16 -> 10 18 eV @. 



III. RESULTS AND DISCUSSION 

Figure shows the total spectrum (Galactic + Extra- 
galactic) compared to the data, assuming continuously 
emitting sources with density n = 10~ 5 Mpc~ 3 and spec- 
tral index 7 = 2.6. The solid and dot-dashed lines for the 
extra-galactic show the median spectrum obtained over 
500 realizations of the source locations. For each real- 
ization the locations of the first hundred closest sources 
(i.e. within ~ 140 Mpc) were drawn at random, using 
a uniform probability law per unit volume; for farther 
sources, the continuous source approximation is valid and 
it was used numerically. The upper and lower dotted 
curves show the 75 th and 25 th percentiles around this 
prediction, meaning that only 25% of spectra are respec- 
tively higher / respectively lower than indicated by these 
curves. This uncertainty is related to the location of the 
closest sources, see below. 

Considering the difficulty of comparing different 
datasets, the fit shown in Fig.n a PP ears satisfying. One 
should also note that this fit uses a minimum number 
of free parameters (7 and Z scat t at 10 17 eV), in order to 
consider the most economical scenario. As discussed be- 
low, there are various ways in which one could extend 
the present analysis, although this comes at the price of 
handling a larger number of (unknown) parameters. 

In Fig. ^ a straight dashed line was drawn across the 
region 1.5 • 10 18 — > 8 • 10 18 eV in which the propagation is 



neither rectilinear nor diffusive. These limits were found 
by comparing the diffusive and rectilinear spectra with 
the no magnetic field spectrum. In this energy range the 
diffusive path length becomes of the same order as the 
rectilinear distance at some point during the particle his- 
tory. The diffusion theorem [14( suggests that the flux 
in this intermediate region should follow the no magnetic 
field spectrum (in which case it would dip ~ 10 % below 
the dashed line around 3- 10 18 eV). This theorem rests on 
the observation that integrating Eq. for a continuous 
distribution of sources over an infinite volume gives the 
rectilinear spectrum Eq. J2J). However the actual volume 
is bounded by the past light cone; this is why the diffusive 
spectrum shuts off exponentially at energies > 10 18 eV. 
The rectilinear part shuts off at energies < 7 • 10 18 eV 
as the maximal lookback time that bounds the integral 
of Eq. J5J) decreases sharply. Hence one might expect a 
small dip in the spectrum around 2 — 3 • 10 18 eV. Interest- 
ingly the data is not inconsistent with such a dip at that 
location. Monte Carlo simulations of particle propaga- 
tion are best suited (and should be performed) to probe 
the spectrum in this region. 

The scattering length was assumed to scale as Z sca tt (X 
E 2 , and its value at 10 18 eV was set here to 17 Mpc. This 
scaling of the scattering length is typical for particles 
with Larmor radius larger than the coherence scale of the 
field, in which case i S catt — Ic(tl/Ic) 2 - Since the Larmor 
radius r L ~ 1 Mpc (E/10 18 eV)(B/10~ 9 G) -1 , one may 
expect this approximation to be valid. In effect, 1 Mpc is 
a strict upper bound to the coherence length of a turbu- 
lent inter-galactic magnetic field [T4I ITEj , and available 
numerical simulations indicate much smaller coherence 
lengths |l6( in clusters of galaxies. A value l c ~ lOkpc 
could also be expected if the inter-galactic magnetic field 
is produced by galactic outflows. The above condition for 
iscatt corresponds to B^/T c ~ 2.5 • 10 -10 G-Mpc 1 / 2 for an 
all-pervading magnetic field. Hence, for l c ~ 20 kpc, and 
B ~ 2 • 10~ 9 G (in order to obtain the correct scattering 
length at 10 18 eV), one finds r L > l c for E > 3 • 10 16 eV. 

It is possible that the scaling of l sca ,tt with energy 
changes in the range 10 16 — > 10 17 eV as may be- 
come smaller than l c . There is no universal scaling for 
^scatt when tl < l c as the exact relationship then de- 
pends on the structure of the magnetic field; for in- 
stance, in Kolmogorov turbulence, one finds i Bca tt oc »"l 

1/3 

for 0.1/ c < tl J$ l c and i Bca tt ^ r i! at lower ener- 
gies 0| • The possible existence of regular components 
of extra-galactic magnetic fields may also modify Z SC att- 
A change in the scaling of Z sca tt with E, if it occurs at 
E > 10 16 ' 5 eV, would imply a different value for B\fQ, 
with the difference being a factor of order unity to a few. 
It is exciting to note that, in the present framework, ex- 
periments such as KASCADE-Grande Q may allow to 
constrain the energy dependence of the scattering length 
(hence the magnetic field structure) by measuring accu- 
rately the energy spectrum and composition between the 
first and second knees. 

The predictions (for both normalizations in Fig. ^) 
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FIG. 2: KASCADE (SYBILL) data: solid lines: all-particle 
spectra (Galactic + extra-galactic); dotted lines: extra- 
galactic component; dashed lines: Galactic spectra, for p, He, 
C, Si and Fe, in order of increasing x— axis intercept. Filled 
circles: all-particle; diamonds: p; upward-pointing triangles: 
He; squares: C; downward pointing triangles: Si; circles: Fe. 
Error bars on reconstucted chemical composition have been 
omitted for clarity, but are significant, see Q. 



for the extra-galactic proton flux are shown and com- 
pared to the chemical composition measurement of KAS- 
CADE in Fig. These composition measurements re- 
main uncertain, as can be seen by comparing the QGS Jet 
and SYBILL reconstructions in 0; the proton and he- 
lium knee positions seem robust however. The dotted 
lines represent the median proton signal from the extra- 
galactic component, whose detection seems within the 
reach of KASCADE-Grande. One may note that Galac- 
tic spectra with exponential suppression beyond the knee 
agree with the KASCADE data. Nonetheless, if the 
Galactic spectra are found to extend as powerlaws be- 
yond the knee, the scattering length of extra-galactic 
protons should be smaller by a factor of order unity (and 
B^fQ correspondingly higher). 

The result for B\fC c depends weakly on the source 
density: since the diffusive (low energy) part of the 
spectrum shuts off as exp [— r 2 /4A 2 ] with r ~ n -1 / 3 
the closest source distance, the cut-off energy de- 
pends on the ratio n _1 / 3 /A oc n~ 1 l z (B\fQ), hence 
Bvt scales with n according to: ByfQ ~ 2 — 3 • 
lO^V/lCT 5 Mpc" 3 ) 1 / 3 G-Mpc 1 / 2 . 

Cosmic variance related to the distance d to the 
closest sources is significant for the low energy (E < 
10 17 eV) and for the high energy (E > 10 20 eV) parts 
of the spectrum, as illustrated by the confidence in- 
tervals around the median flux shown in Fig. ^ hi 
these two energy ranges, the effective linear distance to 
the source is limited to < 50Mpc, which is compara- 
ble to the expected distance to the closest source. In- 
terestingly, the spectra close to both low and high en- 
ergy cut-offs are strongly correlated due to the above 



effect. The distances to within which one should find 
N = (1,2,3,4,5) sources (with N the Poisson aver- 
age) are r ~ (29, 36,41, 46, 50) Mpc respectively. The 
diffusive spectrum sums up contributions that scale as 
exp [— r? /4A 2 ] , with the distance to the i th closest 
source. Therefore, close to the cut-off energy, where 
r 2 /4A 2 3> 1, the total spectrum is dominated on the aver- 
age by the individual spectrum of the closest or the two 
closest sources. At higher energies spectra of more re- 
mote sources contribute with a weight ~ exp [— r 2 /4A 2 ] . 
Since what matters most for the comparison to the data 
is the cut-off energy, one finds that as a first approxima- 
tion, cosmic variance related to the position of the clos- 
est source at distance r*i induces an uncertainty of the 
inferred magnetic field strength AB/B ~ Ad/d ~ 0(1) 
since, as before, the cut-off energy depends on the ratio 
ri/A. 

It is possible that the ultra-high energy cosmic ray 
sources are intermittent with an activity timescale 
^source ^ -^o" 1 ' the previous discussion has assumed 
steady sources corresponding to T sourcc ~ Hq 1 . If 
^source H^ 1 , the number density of sources inferred 
from clustering at high energies (E > 4 • 10 19 eV) un- 
derestimates the actual density of potential sources by 
a factor T/H^ 1 , with T 2 ~ T 2 ourcc + Ar|, and where 
Atb is the typical time spread at E ~ 4 ■ 10 19 eV due 
to magnetic delay. The average time delay reads: tb ~ 
1.5 • 10 8 yr(£;/10 19 eV)- 2 (d/lGpc) 2 (BVZ;/3 ■ lCr 10 G • 
Mpc 1/2 ) 2 0; the ratio At b /t b < 1 depends on the 
structure of the random magnetic field, see Each 
source then contributes for a fraction T/Hq 1 of a Hubble 
time to the diffusive spectrum given in Eq. Jy, but there 
are H^ 1 /T times more sources: the total flux remains 
the same than evaluated previously, except that the cut- 
off energy will correspond to that expected for a source 
density larger by H ~ 1 /T. Hence, following the previ- 
ous discussion, the present scenario remains valid if the 
magnetic field strength is higher by a factor (Hq 1 /T) 1 / 3 . 
For instance, for active galactic nuclei sources of ultra- 
high energy cosmic rays with T source ~ 10 8 yr, a fit 
similar to that shown in Fig. can be obtained for a 
magnetic field strength B\fQ ~ 10~ 9 G ■ Mpc 1 / 2 . For 
the particular case of transient Galactic sources, such as 
7— ray bursts, the situation is different, since the clos- 
est sources lie at distance w 0. Therefore the diffu- 
sive spectrum J dm (E) « £,(c/47r) [47r£>(f?)ti]~ 3/2 q(E), 
since A 2 ~ D(E)ti for close by sources, with ti the look- 
back time to the i th event, and q(E) the injection spec- 
trum per source. Diffusion in extra-magnetic fields thus 
does not produce a low energy cut-off in this case; the 
spectrum is rather subject to the fluctuations of the time 
distribution of past Galactic events. 

At high energies, E > 10 19 eV, particles travel in a 
quasi-rectilinear fashion, i.e. the deflection angle suffered 
by crossing a coherence cell of the magnetic field 69 ~ 
Zc/rL ~ 3 x 10- 3 (Z c /30kpc)(£;/10 19 eV)- 1 (S/10- 9 G) 
is much smaller than unity. The total deflection an- 
gle summed over the trajectory remains smaller than 
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unity, and this justifies the use of Eq. J5J) [l8j : 
6» rms - 25°(£;/10 19 eV)- 1 (d/lGpc) 1 /2(B^/4/3.iO- 10 G- 
Mpc 1 / 2 ). This also implies that charged particle astron- 
omy will be possible at the highest energies. Recent stud- 
ies have attempted to obtain definite predictions for 9 rms 
by using MHD simulations of large-scale structure for- 
mation with magnetic fields scaled to reproduce existing 
data in clusters of galaxies |2Cj . Their results differ 
widely, thereby illustrating the difficulty of constraining 
ah initio the strength of extra-galactic magnetic fields. 
The present value for # rnja _is comparable to or slightly 
larger than that of Ref. [0|j and substantially smaller 
than that of Ref. 20]. The magnitude of 6* rms indi- 
cates that extra-galactic magnetic fields could be probed 
through the angular images of ultra-high energy cosmic 
ray point sources, and this will constitute a strong test 
of the present scenario. 

The proposed scattering length cannot result from 
scattering on magnetic fields associated with galax- 
ies or groups and clusters of galaxies, since the col- 
lision mean free path with either of these objects is 
too large, being ~ O(lGpc). The inferred magnetic 
field might in principle be concentrated around the 
source (on distance scale L) and negligible everywhere 
else. Since the spectrum would cut off below an en- 

1 /2 

ergy such that 2A « 2 [cHq %cattl ^ L, this re- 
quires B > 1/iG (/ c /10kpc)- 1 /2( jL / 10 0kpc)- 1 (for a cut- 
off at 10 17 eV). This possibility cannot be excluded but 
it gives a non-trivial constraint on the source environ- 
ment. Searches for counterparts at the highest energies 
would help test this possibility: for instance, magnetic 
fields such as above are found in clusters of galaxies but 
there is no report of clusters in the arrival directions of 
the highest energy events. If this magnetic field is in- 
trinsic to the source, or if the cut-off at < 10 18 eV is 
due to injection physics in the source y|, then, under 
the present assumptions, the present work still gives a 
stringent upper bound on all-pervading magnetic fields. 
To remain conservative, one may require that the cut- 
off should not occur above ~ 10 18 eV, in which case one 
finds B^/T c < 1(T 9 G-Mpc 1 / 2 . This limit is still an order 
of magnitude below existing Faraday bounds. 

The magnetic field in question thus appears inter- 
galactic in nature, in which case it is likely to be inho- 
mogeneously distributed on small scales. Further studies 
are then required to relate the average B^fQ with the ac- 
tual structure and distribution of these magnetic fields. 
One needs to account for the possible existence of a regu- 
lar magnetic field component aligned with filaments and 
walls, which would inhibit perpendicular transport, and 
consider the respective filling fractions and amplitudes of 
the turbulent and regular components. It would be cer- 
tainly worthwhile to extend the simulations of particle 
propagation made in realistic magnetic fields [1(1 [20] to 
the energies of interest. 

Finally there are various ways in which the present 
study could be extended. One should notably consider 



the possible energy dependences of the scattering length 
(including the above effects of inhomogeneous and regu- 
lar magnetic fields), the role of intermittent sources, the 
possible cosmological evolution of the magnetic field and 
of the source density, and, as mentioned above, the pos- 
sibly inhomogeneous structure of the magnetic field on a 
scale comparable to the closest ultra-high energy cosmic 
ray sources. 
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APPENDIX: DIFFUSION OVER 
COSMOLOGICAL SCALES 

The diffusion of particles in an expanding background 
space-time can be seen as a standard diffusion process on 
a fixed background in conformal coordinates (ry, r) , with 
r] the conformal time defined by a(ri)dri = dt and r the co- 
moving coordinates in a Friedman-Lemaitre-Robertson- 
Walker metric; a denotes the scale factor and t cosmic 
time. One can indeed approximate the diffusing process 
as a random walk against scattering centers of constant 
comoving coordinates. 

Particles also experience dilution due to expansion, 
expansion energy losses and energy losses due to pair 
and pion production on diffuse backgrounds. At redshift 
z = 0, these photo-interaction losses are negligible with 
respect to expansion losses for energies E < 2 ■ 10 18 eV, 
but become increasingly more important at higher red- 
shift due to the increased cosmic microwave background 
temperature and density Nonetheless the main en- 
ergy loss in the course of the history of a particle with 
present energy Eo < 10 18 eV is due to expansion. One 
reason is that the majority of the sources that contribute 
to the diffuse flux at energy Eg are located at moder- 
ate redshifts as a result of the nonlinear time-redshift 
relation: redshift z = 2, for instance, corresponds to a 
lookback time of 76% of the age of the Universe [2l| . 
More importantly, pion and pair production losses at 
high redshift become catastrophic, so that the time in- 
terval during which the losses are dominated by photo- 
interactions is much smaller than a Hubble time. Fi- 
nally, the contribution to the diffuse flux at energy Eq 
of particles injected with energy E g scales, in a first ap- 
proximation, as q(E g )dE g /dE, with q(E s ) the injection 
spectrum and dE g /dE acounts for the dilation of the en- 
ergy interval. The function E s (rj,E) defines the energy 
of the particle at time ?y, assuming it has energy E at 
time to- This function and its derivative dE g /dE can be 
reconstructed by integrating the energy losses 0. As- 
suming dE g /dE « E g / E, which is exact for expansion 
losses, one sees that the contribution of particles injected 
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at remote lookback times (hence with high E g ) is negligi- 
ble with respect to that of particles injected recently with 
E g « E since q(E g ) cx E~ J and 7 ~ 2.6 here. The nu- 
merical difference between a diffuse flux computed using 
only expansion losses and that computed with all energy 
losses included is indeed less than 5% at E -C 10 18 eV, 
and increases to 20% at E ~ 10 18 cV. Consequently, it 
is assumed in this discussion that particles with present 
-Eo 10 18 eV have been subject to expansion losses only 
throughout their history. 

Energy losses due to expansion are expressed as: 
dE/dr] = —HE, with Ti = (l/a)da/drj the expansion 
rate in conformal time. The phase space density of 
particles N(r], E, r) at coordinates r, time 77 and en- 
ergy E, which is related to the distribution function by 
N(t},E,t) — (Anp 2 / c) f (r] , p, r), with pc — E, is solution 
to the diffusion equation: 



d_ 

drj 



d 



a 3 N)-VDV (a 3 N)-H— (Ea 3 N) = a 3 Q (ry, E, r) . 

(A.l) 

where Q gives the (physical) number density of particles 
injected per unit energy and conformal time intervals. 
The a 3 prefactor of N takes into account the effect of 
dilution of particle density through expansion. The fact 



that expansion losses are separable in terms of the two 
variables E and rj allows to find an exact solution to 
this diffusion equation, using standard Green functions 
methods; see |22( for the same problem with time inde- 
pendent losses in a non-expanding background. Explic- 
itly, through the change of variables (77, E) — > (u, v), with 
u = log(aE) and v = log(a/E), one can derive the Green 
function (for the equation for N) as: 



G (r) , E Q ,r ;r) e , E e ,r e ) 



exp 



Iro-r.l 2 



(4ttA 2 ) 



3/2 



* f E. _ ^ I , (A.2) 



with the shorthand notations: do 
The path length A is defined by: 



0(770) and a e = a(f7 e 



A 2 



no 



dr/D 



CLpEp 



a(rj) 



(A.3) 



where D(E) is the diffusion coefficient; if this latter de- 
pends explicitly on time, for instance if the magnetic field 
strength evolves with redshift, the solution remains valid. 
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